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ABSTRACT 

Type la supernovae (SNIa) remain mysterious despite their central importance in cosmology and their rapidly 
increasing discovery rate. The progenitors of SNIa can be probed by the delay time between progenitor birth 
and explosion as SNIa. The explosions and progenitors of SNIa can be probed by MeV nuclear gamma rays 
emitted in the decays of radioactive nickel and cobalt into iron. We compare the cosmic star formation and 
SNIa rates, finding that their different redshift evolution requires a large fraction of SNIa to have large delay 
times. A delay time distribution of the form t^" with a = 1.0 ± 0.3 provides a good fit, implying 50% of 
SNIa explode more than ~ 1 Gyr after progenitor birth. The extrapolation of the cosmic SNIa rate to z = 
agrees with the rate we deduce from catalogs of local SNIa. We investigate prospects for gamma-ray telescopes 
to exploit the facts that escaping gamma rays directly reveal the power source of SNIa and uniquely provide 
tomography of the expanding ejecta. We find large improvements relative to earlier studies by Gehrels et al. in 
1987 and Timmes & Woosley in 1997 due to larger and more certain SNIa rates and advances in gamma-ray 
detectors. The proposed Advanced Compton Telescope, with a narrow-line sensitivity ^ 60 times better than 
that of current satellites, would, on an annual basis, detect up to ^ 100 SNIa (Scr) and provide revolutionary 
model discrimination for SNIa within 20 Mpc, with gamma-ray light curves measured with ^ 10a significance 
daily for ^ 100 days. Even more modest improvements in detector sensitivity would open a new and invaluable 
astronomy with frequent SNIa gamma-ray detections. 

Subject headings: Gamma rays: diffuse background — Gamma rays: stars — Nuclear reactions, nucleosyn- 
thesis, abundances — Stars: statistics — Supernovae: general 



1. INTRODUCTION 

Type la supernovae (SNIa) are deeply connected with 
many important frontiers of astrophysics and cosmology. 
They occur in all galaxy types and are major contributors 
to galactic chemical evolu tion, in particular of iron (e.g., 
iMatteucci & Greggiolll98d) . They are very bright and, as 
high redshift distance indicators, play a cri tical ro le in estab- 
lishing the modern cosmology paradigm (Riess et al. 1998; 
iPerlmutter et"ai]|l999l:lFreedman et ali200liK 

However, there are major uncertainties regarding the nature 
of the SNIa progenitors and explosions. While it is estab- 
lished that most SNIa result from the thermonuclear explosion 
of carbon-oxygen white dwarfs (WD) near the Chandrasekhar 
mass, the mechanism of mass gain remains debated. In the 
single-degenerate (SD) sc enario the WD accretes mass from 
a companion star dWhelan & Iben 1973) . while in the double- 
degenerate (DD) scena rio the WD merges with another WD 
(llben&Tutukovlll98 4'). In addition, although the main prod- 
ucts are known, the basic mechanism of nuclear burning re- 
mains under debate. In deflagration the ignited flame prop- 
agates subsonically, while in detonation it propagates super- 
sonically as a shock wave. There are combined models, as 
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well as the possibility of deto nation in the He layer of a sub- 
Chandrasekhar WD (see, e .g., lHillebrandt & Niemeveiil2000l : 
iTutukov & Fedoroval200"7l) . 

The difference between star formation and SNIa rates de- 
pends on what the SNIa progenitors are. The rate of SNIa fol- 
lowing an episode of star formation depends on (i) the delay 
time, which describes the time required for a newly-formed 
SNIa progenitor to develop into a SNIa, and (ii) the produc- 
tion efficiency, which is the number of SNIa produced per 
Mq of star formed (e.g., Greggio 2005). The SD and DD 
scenarios involve different progenitors and have different pre- 
dictions for these quantities. Although current predictions 
are comparable within uncertainties, improved information 
on the global delay-time distribution (DTD) and SNIa effi- 
ciency will discriminate b etween SNIa progenitor scenarios 
dYungelson & Livioll2000l ). 

It is well known from nuclear fusion modeling and opti- 
cal light curve observations that each SNIa yields as much 
as 0.5-0.7 Mq of radioactive ^"^Ni. The decay of ^""Ni, via 
^^Co to stable ^^Fe, provides the primary source of energy — 
in the form of nuclear gamma rays and energetic positrons that 
deposit much of their energy in the ejecta — that powers the 
SNIa optical display. Initially, the gamma-rays are trapped, 
but as the SNIa ejecta expands and the matter density drops 
they start to escape, while the positrons remain largely trapped 
until much later times. The decline rates of the optical light 
curve show the timescales set by the decay of ^^Ni (half life 
of 6 days) and '"^Co (half Ufe of 77 days). 

The detection of gamma rays that escape th e ejecta is the 
key to resolving the central mysteries of SNIa jGehrels et alj 
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Il994t iHoflichet al.lll998t lBoggsll2006l) . Since gamma rays 
are more penetrating and their opacities are much simpler 
than those for optical photons, they offer a more straight- 
forward and direct probe of the inner mechanisms of SNIa 
jClavton et al.ll969HMilne et alj2004irsern et alj2008l) . The 
gamma ray flux allows the identification of radioactive mate- 
rial yield, and its time evolution allows tomography of the 
surrounding ejecta; both quantities can be compared with 
model predictions. In a few cases of fortuitously nearby 
SNIa, hmits just above theoretical expectations have been se t 
(lMatz&Sharelll990l: iLichti et al.lfT99i iLeising et alJll999l) . 
No SNIa features have been seen in the co smic gamma-ray 
background (CGB; e.g. JStrigari et al.ll2005h . ReHable knowl- 
edge of SNIa rates is essential to defining realistic prospects 
for gamma-ray detection, and new data make this possible . 

In the first part of this paper, we investigate the cosmic SNIa 
rate and determine the delay of SNIa relative to their progen- 
itor formation and the efficiency of forming SNIa (Section 
|2|. We use a large selection of star formation indicators with 
recent updates (Section lzTl ). as well as a comprehensive com- 
pilation of SNIa rate data (Section [2.2l l. We study the effects 
of delay and efficiency using data over a substantial redshift 
range, and we discuss what our results reveal about the SNIa 
progenitors, i.e., for SD and DD scenarios (Section l23] l. 

In the second part of this paper, we bring in gamma rays as 
a probe of SNIa explosions and progenitors (Section|3]l. Our 
analysis of the cosmic SNIa rate plays a critical role in deter- 
mining the gamma-ray detection prospects. We first review 
the gamma-ray emission for several benchmark SNIa models 
(Section lTTl l. We determine the SNIa contribution to the CGB 
using our updated cosmic SNIa rate inputs (Section [J!2T l. We 
derive SNIa rates from supernova catalogs, and joining them 
with cosmic SNIa rates, we investigate the local SNIa rate 
with a higher accuracy than previously possible (Section [373] i. 
We revisit the status of gamma-ray observations in light of 
our updated local SNIa rate (Section 13.4b . Finally, we give 
new results on the excellent gamma-ray detection prospects 
for nearby SNIa and the physics that this will probe (Section 

We close with a summary of our findings (Section |4|l. 
Throughout, we adopt the standard ACDM cosmology with 
n,n = 0.3, = 0.7, and Hq = 70 km s^^ Mpc^^. 

2. SNIA PROGENITORS AND COSMIC SNIA RATE 

Comparing the cosmic star formation rate and the high red- 
shift SNIa rate of the GOO DS fields sug gests a DTD tightly 
distributed around 3-4 Gyr dStrolger et a l. 2004; Dahlen et akl 
l2008t IStrolger et al.ll2010h . However, small SNIa statistics, 
uncertain dust corrections, and the uncertain star formation 
rate at high redshift comp licate studies of the DTD in this 
range (Greggio et al. 2008). We perform a new investigation 
of the SNIa rate by adopting both greatly updated star forma- 
tion rate data (Section [2. lb and SNIa rate data (Section [2.2b . 
In particular, we focus on the low redshift range, which pro- 
vides a complementary test to previous studies. By comparing 
the calculation to observations, we determine the global DTD 
and SNIa efficiency, and discuss implications for SNIa pro- 
genitors (Section l23T i. 
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Fig. 1. — Comoving cosmic star formation rate density as a function of 
(1 + z) on a logarithmic axis. Data as labeled (see text for details), agree 
to within ±20% for z <\. The uncertainty increases with redshift, but the 
overall shape is we ll constrained. The fitted curve is shown by a solid line 
iYuksel et al.l2008l) . 

2.1. SNIa rate calculation 

The SNIa rate at a given time corresponds to the stellar 
birthrate at earlier epochs, convolved with the DTD, normal- 
ized by the efficiency of forming SNIa progenitor systems. 
The normalization is 



(1) 



?7= / Aia(Af)e(M)dM, 



where Aia(Af) is the SNIa formation efficiency, ^(Af ) is the 
initial mass function, and the integration is performed over the 
mass range for SNIa progenitors. Unless otherwise stated, we 
present all results for the Salpeter initial mass function; our 
results on the SNIa rate do not strongly depend on this choice 
because the SNIa progenitor mass range is sufficiently high. 
Assuming Aja and ^ do not vary with space and time, 77 is a 
constant representing the number of SNIa produced per Mq 
of stars formed, so that 



Jtxa 



(2) 



where is the DTD, its argument is the time difference be- 
tween progenitor star formation and SNIa explosion, p*(tc) 
is the star formation rate, and the integration is performed 
over cosmic time from tio, the age of the universe when 
the first stars were being formed, z ~ 10. is normal- 

ized such that its integral over the entire delay time range 
is unity. We do not con sider the effect of metallicity on 
the rates (e.g., Prantzos & Boissier 20 03; Prieto et al, 2008[ 
iBoissier & Prantzos[|2009l;lMeng et al.ll2009i) . 

In the limit that $(<) is large only for small t, we recover a 
prompt population of supernovae, i.e., similar to core-collapse 
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supernovae, for which the cosmic rate evolution is the same 
as for the star formation rate. For SNIa, the observation that 
they occur in both young star-forming galaxies and old ellip- 
tical galaxies with little star formation demonstrates the need 
for a wide range of delay times. Recently, two populations 
of SNIa, a "prompt" population with delay times ~ 0.1 Gyr 
and a "delaye d" population w ith delay ti mes ^10 Gyr, have 



been r eported (IMannucci et al.i 2005; Sc annapieco & BildstenI 
20051: iMannucci et al.ll2006t ISulhvan et al.ll2006l: iMaoz et al.l 



2010ah . Whether this represents a truly bimodal distribution, 
or whether it is the result of two coarsely sampled bins of what 
is in reality a continuous DTD, remains uncertain. 

We adopt a power-law DTD, cx t^". Power-laws have 
been reported from SNIa observation s: for exa mple, with 
a = 0.5 ± 0.2 from the SNLS data (Pritchet et al.li2008.) . with 
a = l.OS^Q fr om observations of transients in passively 
evolving galaxies ('Totani et al.l 1200 8*), with a ~ 1 from lo- 
cal SNIa in the Lic k Observatory Supernova Search (LOSS) 
dMaoz et al.ll2010ah. and most r ecently a = 1.2 ± 0.3 from 
cluster SNIa ("iVIaoz et al.ll2010bl) . In fact, the number of WDs 
following an instantaneous star formation burst is also a sim- 
ple power-law in t: for an initial mass function £,{M) oc 
(x = —2.35 for Salpeter) and a power-law approximation 
to the evolutionary time-scale of 3-8 M© stars tagc « 
(y w —2.6), the index isQ: = {x — y + l)/y ^ 0.5. However, 
this distribution extends only over the lifetimes of progenitors, 
30^00 Myr, which is much smaller than the observed range 
of delay times. The extra time needed for mass transfer be- 
fore triggering a SNIa would shift the delay time distribution 
towards longer delays, but the final DTD is largely dependent 
on the properties of the binary and its evolution. 

We apply our DTD over times 0. 1-13 Gyr and adopt a as a 
parameter. The lower time limit could be as short as the Ufe- 
time of a 8 Mq star, ~ 0.03 Gyr, but recent studies of the spa- 
tial distribution of SNIa in their hos t galaxies indicate a lower 
Umit of 0.2 Gyr (iRaskin et al.l2009l) . We adopt 0.1 Gyr; at the 
precision of the current SNIa rate measurements, our results 
are not very sensitive to small delays. We fix our upper limit 
to the look-back time to z = 10, which is sufficiently long 
to accommodate SNIa observed in the oldest elliptical galax- 
ies. Allowing a larger value would not change the redshift 
evolution, and would simply scale the efficiency. 

The cosmic star formation r ate density is shown in Fi: 
ure [T] The compilation of iHopkins & BeacomI (1200^ 
which includes data from multiple star formation indi- 
cators, is shown togethe r with recent data fro m mid- 
infrared-emitting galaxies ( Ruiopakarn et al.l 2010l). Lyman 



break galaxies (LBG) ( Verma et alj l2007l iBouwens et al.l 
1^07, 2008; Reddy & Steidel 20oC and gamma-ray burst 
(GRB) host galaxies dVuksel et al. 20081: iKistler et all l2009l 
Wanderman & PiranI I20r0t iButler et alTbOl d). Note that 



Bouwens et al. (2007 2008b report results integrated down 



to a chosen lu minosity limit ; we sh ow instead complete in- 
tegrations (see iKistler et al.l (120091) for details). We also 
show Hubble U ltra Deep field (UDF) measurements from 
lYan et al. 1 (1201 Oi) . similarly integrated to low galaxy lumi- 
nosities. New UDF results from Bouwens et al.l (l20ld) have 
just appeared. As these indicate a steeper rise at low galaxy 



luminosities, the integration for the full star formation rate 
would need to be carefully regulated, and we do not perform 
this here. The star formation rate has been checked by in- 
dependently meas ured quantities, such as the extragalactic 
background light (Horiuchi et al.ll2009h . stellar mass density 
dWilkins et al.il2008h . and upper limits on the diffuse super- 
nova neutrino background (Beacom 2010). 

On the top axis we label cosmic time. It highlights that 
delay times of Gyrs will result in noticeable differences in the 
slope between the star formation and SNIa rate rates. The star 
formation rate rapidly increases by one order of magnitude 
between redshift and 1, which has been very well measured 
by a variety of indicators. This serves as an important feature 
to be compared to the SNIa rate. At higher redshift, the star 
formation rate eventually declines, at first slowly after z ~ I, 
and then more rapidly after z ^ 4. The exact slope above 
z ^ A remains uncertain, although this does not affect our 
conclusions, as we discuss in Section l23] 



2.2. SNIa rate measurements 

Since the first measur ement of SNIa at cosmological dis- 
tances dPainet all [19961) . many surveys have collected ho- 
mogenous samples of CCD-discovered SNIa. These surveys 
periodically observe the same patch of sky, or the same sam- 
ple of galaxies, to find transients. Cuts are applied to select 
the most confidently-identified supernovae, and corrections 
for dust and incompleteness are applied. We summarize the 
present SNIa rate measurements in Table |3] 

The high redshift {z > 0.5) rates are the hardest to mea- 
sure. Dust corrections are non-trivi al and uncertain, e.g., the 
dust corrected data of iDahlen et al. I f2004l2008h are a factor 
^ 2 higher than their uncorrected data. Even this has been ar- 
gued t o be insufficient and the rates could be treated as lower 
Umits dGreggioetal .112008 ). On the other hand, SNIa sur- 
veys are likely contaminated by core-collapse supernovae, for 
which the rate increases from « 10^"* yi'^^ Mpc^'^ at z = 
to w 10"^ yr^i Mpc"^ by z = 1, i.e., 10 times the SNIa 
rate by z = 1. It is therefore important that only the most 
confidently-identified SNIa — those that are spectroscopically 
confirmed — are selected for rate measurements. Indeed, 
rates derived using photome t ric classification ( Bjmig & Tonr^ 
2006; iPoznanski et al.ll2007l iKuznetsova et al.ll2008l) are gen- 
erally higher than those bas ed on a h igher fraction of spectro- 
scopically confirmed SNIa dDahlen et al. 2004, 2008). This is 
usually reflected in the large systematic errors. Finally, the 
sample sizes rapidly become small at high redshift, where 
typically A^ia = 2-3. These issues make it difficult to use 
the high redshift rate measurements dOda et al.ll2008b . even 
though they would be very sensitive to delays due to the small 
cosmic time difference since the first stars were being formed. 

At intermediate redshift (0.1 ^ z < 0.5), many of the 
measu rements are based solely on spectroscopic confirmed 
SNIa (iHardin et al.ll2000l: iPain et al. II2002I: iTonrv etani2003l: 



Madgwick et al.' '2003t iBlanc et al.1 120041; iNeill et al.1 12OO6I; 
Dilday et al. 2008). Multi-band phot ometric identification is 
also now rehable: iDildav et al.1 d2010l) show that in the SDSS- 
II Supernova Survey only ^2% of the photometric SNIa in 
z < 0.3 may be misidentified. Furthermore, the effects of dust 
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Fig. 2. — Comoving SNIa rate as a function of redshift. Measurements are 
shown with statistical and systematic uncertainties combined in quadrature. 
The filled po ints are more reliable than the empty points (see text). The local 
ISmartt et alj J2009j) point is affected by a local enhancement (see Section 
13. 3t . The best-fit a = 1.0 power-law DTD (solid) is shown, with a shaded 
band reflecting the uncertainty. Other DTD shown are as labeled in Figure[3l 
they do not fit the data as well. 

become less important than at high redshifts dBotticella et alJ 
I2OO8,) . However, many surveys target a pre-selected sample 
of galaxies and, while a large sample is adopted, there would 
be a bias towards bright galaxies, and supernovae that explode 
in faint galaxies would be missed. Lastly, we note that sample 
sizes vary largely from survey to survey. 

Finally, in the local {D < 100 Mpc) volume, the rate 
has be en calculated from compilations of carefu lly selected 
SNla dCappellaro et al ] |l999t ISmartt et al. 2009). as well as 
the local LOSS measurement (Leaman et al., ,2010t iLi et aT 



l2010allbh . Note that the rate calculated from ISmartt et al" 
( l2009h probably reflects a local enhancement (Section [33T l 



2.3. SNIa rate synthesis 

Compared to the star formation rate, the SNIa rate is fur- 
ther away from a consensus. Although the number of rate 
measurements has increased, the scatter is large, and at times, 
measurements in comparable redshifts are in direct disagree- 
ment. Bearing in mind the strengths and weaknesses of the 
data, we group the rate measurements into two categories: 

• Filled points: those based on dedicated surveys using 
only spectroscopically identified SNIa. We also in- 
clude select surveys with less than 100% spectroscopic 
identification, including the survey by SDSS which 
includes an order of magnitude more SNIa than any 
other survey (Dilday et al. 2010), and measurements by 
lOahlen et all (|2004 l2008l) which are the most reliable 
measurements in the high redshift regime. 

• Empty points: those based largely on photometrically 
identified SNIa. 
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Fig. 3. — Same as Figure m but over the entire redshift range where data 
are available. In addition to the best-fit a = 1.0 power-law DTD (solid), 
we show the no delay (dotted), the two-component DTD of Mannucci et aU 
flood) (dot-dashed), and the 3.4 Gyr naiTow DTD of iDahlen et al., t2Q08l) 
(dot-dot-dashed), as labeled. 



We caution that even within each category, the sample size, 
observation schedule, limiting magnitude, and other condi- 
tions vary, and our description is simply an attempt to appreci- 
ate so me of t he important differences between measurements 
jBlanc & Greggioll2008h . 

Now, we can make better sense of the calculated and mea- 
sured SNIa rates despite their uncertainties. When we restrict 
ourselves to the filled points, we find that the SNIa rate in- 
creases by only a factor ^ 2-3 from redshift to 0.8 or so. 
On the other hand, the cosmic star formation rate increases 
by a factor 10 from redshift to 1. This implies that a large 
fraction of SNIa have cosmologically-large time delays. 

We show in Figures |2] and [3] the calculated and measured 
SNIa rates. We fit the calculated cosmic SNIa rate form to 
the selection of rate data consistin g of filled po i nts an d local 
measurements (excluding that by Smartt etaTI (l2009h ). For 
the power-law DTD we fit for the normalization at redshift 
zero, i?ia(0), and the DTD exponent, a. We find best-fit val- 
ues (i?ia(0),a) = (0.24, 1.0) with = 5.4 for 17 degrees of 
freedom. The projections of the Ax^ 2.30 elliptical con- 
tour give i?ia(0) = 0.24 ± 0.04 and a = 1.0 ± 0.3; these 
parameters are strongly anti-correlated. The uncertainty is 
shown as grey shading in the figures. The best-fit efficiency 



and uncertainty of making SNIa is (5 ± 1) x 10~ Mq , or, 
assuming a SNIa progenitor mass range of 3-8 M0, a SNIa 
fraction (which we define as the number of SNIa divided by 
the number of 3-8 Mq stars formed; the fraction of stars in 
SNIa-producing binaries is twice this) of 2.4 ± 0.5%. 

In add ition to the pow e r-law DTD we show the bimodal 
DTD of iMannucci et all (l2006h (dot - dashed), the 3 4 Gyr 
naiTow DTD of 'Str olger et al.1 (I2OO4I) ; iDahlen et all (l2008b ; 
iStrolger et al.. (.2010) (dot-dot-dashed), and the no-delay case 
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(dotted), for comparison. These DTDs do not fit the data as 
well, with values of 10, 17, and 21, respectively. All rise 
too fast compared to the reliable data (filled points). This also 
m eans they unde r predic t the z = rate. The d ifferent results 
in IStrolger et all (12004 : iDahlen et all (l2008h : IStrolger et all 
(1201 Ol) are due to the choice of data. The narrow DTD is 
driven by a declining SNla rate at z ^ 1.4 (see green left- 
pointing triangles in Figure |3]l, while our result is driven by 
the slow rise at low redshift. 

At present, the DTD does not clearly identify the SNla pro- 
genitor scenario, because both the SD and DD scenarios can 
lead to power-law DTDs with a w 1. In the DD scenario, the 
delay is approximately dictated by the time taken by a binary 
to merge by angular momentum loss, which from general rel- 
ativity scales as the fourth power of the binary separation. For 
a logarithmically-flat distribution of binary separations, one 
obtains a ^ 1. Although binary evolution synthesis stud- 
ies of the SD scenario typically show that the DTD peaks 
at characteristic scales related to stable rnass transfer (e.g., 
Yungelson & Livio"200(f; 'Bel czvnski etaLll2005l iMeng et all 
2009; Meng & Yang 2010 ), power- law DTDs have also been 
predicted: for example, iHachisu et al.l ( l2008l) report power- 
law DTDs with awl. Both progenitor scenarios are ac- 
ceptable fits, but as predictions of SNla progenitor scenarios 
improve and SNla rate data accumulates, better testing would 
become possible. We have assumed for simplicity a power- 
law DTD: in order to generally test the progenitor scenario, a 
generic form for the DTD should be statistically tested (e.g., 
IStrolgeretalhoiOl) . 

The choice of the initial mass function does not affect 
the shape of the SNla rate, but changes the number of 
SNla progenitor stars. This introduces a small tens of per- 
cent difference in the efficiency. For a modern Baldry- 
Glaz ebrook init ial mass function with a low-mass suppres- 
sion dBaldrv & G lazebrook 2003), the SNla fraction required 
is 2.9 ± 0.6%. Since the comparison of cosmic rates is not 
sensitive to small delays, there could be more short-lived pro- 
genitors than indicated by our simple DTD form. Uncertainty 
in the cosmic star formation rate beyond z ~ 4 affects the low 
redshift SNla slope, albeit weakly, through the long tail of the 
DTD. There is a degeneracy in the high-redshift star forma- 
tion slope and the value of a, where a higher star formation 
rate can be compensated by a larger a. However, at the preci- 
sion of the current SNla rate measurements, this uncertainty 
does not affect our conclusions. 

Our SNla efficiency is compara ble to those from recen t 
studies of volumetric SNla rates dBlanc & Greggiol l2008h . 
and somewhat lower than those derived by other methods 
lOreggio 2010: Maoz & Badenes 2010: Maoz et al. 2010a|). 
Note that comparisons must be done with a common assump- 
tion for the initial mass function. As described above, a more 
modern initial mass function would increase our efficiency, 
although a difference remains. The true value o f the effi- 
ciency remains to be determined: see ^ Maozl (l2008i) for a de- 
tailed discussion and implications of the SNla efficiency. For 
example, our SNla fraction contrasts with the ^ 0.2% advo- 
cated in some DD scenarios, where conservativ e conditions 
on the mass ratios are applied dRuiter et al.ll2009h . The higher 



required efficiency that we found suggests the need to relax 
some assumptions, or to include contributions from SNla of 
lower mass progenitors. We should note here that our def- 
inition for efficiency implicitly assumes binaries that would 
become a SNla within w 13 Gyr. 

It is revealing to discuss the SNla rate jointly with the core- 
collapse supernova rate. The ratio of SNla to core-collapse 
supernova (la/11) at 2; = is « 0.2-0.3. Here we use the 
z ~ cosmic SNla rate instead of the local 10 Mpc data, 
as the latter rate is ~ 0.1 per year: the r ecent decade, for 
example, had no SNla dKistler et al.ll2008b . Since SNla are 
delayed, the SNla rate should increase less strongly than the 
star formation rate and the core-collapse rate from z = to 
z — 1. The la/11 ratio should therefore decrease with redshift 
in this range. However, the data do not cle arly show such 
a tren d, being consistent with no evolution (IHoriuchi et al.l 
("200^: this is true even with updates on the core-collapse rate, 
e.g., .Bazin et a 1] (l2009b ). A Ukely possibility is that fainter 
core-collapse supe rnovae are being mis sed compared to the 
brighter SNla (see iHoriuchi et alll2009l) . The importance of 
the missing core-collapse supernovae spans topics such as the 
formation of black holes in massive stellar collapse, super- 
nova neutrinos, and metal enrichment. 

3. SNIA EXPLOSIONS AND GAMMA-RAYS 

Gamma rays reveal information of the SNla interior that can 
be effectively used to study SNla physics. We first review the 
gamma-ray emission per SNla, focusing on a range of SNla 
models that accommodate normal, superluminous, and sublu- 
minous SNla (Section [3. It . Using the cosmic SNla rate an- 
alyzed in Section 12] we discuss prospects for using the CGB 
to detect SNla gamma rays (Section [J!2T i. We then discuss the 
local SNla rate with higher precision than previously possible, 
by combining our cosmic SNla rate analysis with SNla rates 
from supernova catalogs of nearby supernovae (Section [3.3t . 
After reviewing the current status of gamma-ray observations 
(Section [3 .4b . we give new results detailing the prospects for 
SNla gamma-ray detection and studying SNla physics. 

3.1. Gamma-ray yield per SNla 

In ah SNla models, the decay chain 56Ni ^ ^^Co^ ^^Fe 
provides the primary source of energy that powers the SNla 
optical display. The ^^Ni decays by electron capture and the 
daughter ^^Co emits gamma rays by the nuclear de-excitation 
process 

^^m + e- ^^^Co* + Ue (3) 

= 158keV(99%), 812keV(86%), 

where percentages express photons per decay and the sum can 
be larger than 100%. Only the dominant lines are noted here. 
The daughter ^^Co decays by electron capture (8 1 %) as well 
as positron emission (19%), and the daughter ^^Fe de-excites 
by emitting gamma rays 

^^Co + e-^'^^Fe* + i^e (4) 
^^Co^^'^Fe* +e+ +i^e 
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TABLE 1 

Gamma-ray line emission from SNIa 











Line flux over 10" s [10"" s"^] 


SNIa model 




Ref. 


AfNi [Mq] 


812 keV 


847 keV 1238 keV 


Normal: delayed-detonation 


DD202C 


1 


0.72 


0.8 


5.2 3.9 


Normal: deflagration 


W7 


2 


0.58 


0.5 


4.3 3.2 


Normal: He-detonation 


HED8 


3 


0.51 


1.6 


4.4 


Superluminous : late-detonation 


W7DT 


4 


0.76 


2.2 


5.9 


Superluminous: He-detonation 


HECD 


5 


0.72 


2.5 


5.5 


Subluminous: He-detonation 


HED6 


3 


0.26 


0.6 


2.2 


Subluminous: pulsed delayed-detonation 


PDD54 


6 


0.14 


0.05 


1.2 
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= 847keV(100%), 1238 keV(67%), 

where the percentages include the effects of the 19% branch- 
ing ratio for positron production by beta decay. 

Initially, most of the gamma rays Compton scatter on the 
SNIa ejecta and deposit their energy. As the ejecta expands 
and the matter density decreases, more gamma rays escape. 
The observed 847 keV line luminosity is 



S^{t) = pcsc(i) 



M^Na 



56 



-t/rca 



't/TNi 



TCo — TNi 



7S 



(5) 

where pcsc is the model-dependent probability of escape 
through the SNIa ejecta, the quantity in the first square brack- 
ets defines the model-dependent nucleosynthesis yield, and 
the second square brackets reflect the nuclear decay rate. 
Here, A/ni is the nickel mass, Na is the Avogadro number, 
and TNi — ti/2/ln(2) where ti/2 is the half life: ti/2 = 
6.1 days for ^^Ni and 77.2 days for '"^^Co. The distribution 
ofM^i ( over many SNIa ) and especially the time- evolution of 
Peso (per SNIa) can be used to probe SNIa physics. 

In deflagration, nuclear burning ignites near the center and 
burning moves subsonically across the progenitor In pure 
detonation, the flame front propagates supersonically as a 
shock front, but we do not consider this further s ince the re- 
sultin g elemental abundances disagree with data CArnett et al.l 
Il97lh . In delayed-detonation, an initial deflagration becomes 
a detonation at some critical density that is an unknown pa- 
rameter. All these models are usually assumed to be initi- 
ated from a WD near the Chandrasekhar mass, accreting mass 
from a non-degenerate star. In He-detonation, burning com- 
mences in the degenerate helium layer near the surface of a 
sub-Chandrasekhar mass WD. In DD scenario models, igni- 
tion occurs at low densities as mass from the disrupted binary 
accretes. It is expected that a large material envelope covers 
the burning sites, so that the escape probability for gamma 
rays is lower than in SD scenarios, thus probing the progeni- 
tors; 

iMilne etakl (12004*) compared the gamma-ray emission from 
seven transport codes for a selection of SNIa models. They 
conclude that differences due to transport codes are 10-20%, 
much less than the differences that result from models. The 
continuum emission is more model dependent than the line 
emission since it depends on multiple scatterings and the 
time-integrated continuum differs by up to a factor ^ 2 be- 
tween codes. 



From lMilne et al J (12004') we adopt a selection of SNIa mod- 
els representative of normal, superluminous, and subluminous 
SNIa, and spanning the deflagration, delayed-detonation, and 
He-detonation models. The average peak line emission, over 
a 10^ s period, are summarized in Table [T] The light curves 
for super-luminous SNIa (solid) and normal SNIa (dashed) 
are shown in Figure [T] We comment on DD SNIa models 
in Section 13.61 We adopt two time-i ntegrated gamma-ra y 
spectra: a deflagration model W7 of iNomoto et alj (Il984l) . 
which yields which yields 0.58 M(7) of ^^Ni, a nd a d elayed- 
detonation model 5p0z22.23 of iHoflich etal] (|2002|) . which 
yields O.56M0 of ^^Ni. Both models are representative of 
normal SNIa, the most common kind. 

3.2. Contribution to the cosmic gamma-ray background 

The SNIa contribution to the CGB depends on the cosmic 
SNIa rate and the time-integrated gamma-ray number spec- 
trum per SNIa, f{E), as 



E'- 



-J{E') 



(6) 



(1 + ^) 

where E is the measured photon energy, and |dz/dt| = 
Hq{1 + z)[an(l + zf + ^Kf/"^. The left-hand side is equiv- 
alent to i^/y, and the redshift factor in Eq. ^ comes from the 
energy scaling. 
100 
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Fig. 4. — The SNIa contribution to the cosmic gamma-ray background. 
The contribution from the deflagration model W7 (solid) and the delayed- 
detonation model 5p0z22.23 (dot-dashed) are shown. For W7 we show the 
range owing to the uncertainty range of the SNIa rate (shaded); a similar 
range applies for 5p0z22.23 but is not shown. Data are as labeled. 
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Fig. 5. — Number of SNIa in the SAl catalog in 10 year bins, plotted differ- 
entially in distance, with selected SNIa as labeled. The solid step shows the 
cosmic SNIa rate extrapolated to the local volume, revealing the incomplete- 
ness of the catalog at large distances. On the top x-axis, the distance is con- 
verted to a 847 keV line flux, assuming 4.3 X 10'*'^ s^^ at peak. The dashed 
lines represent Sa detector line sensitivities (10*^ s), labeled with years of 
operation; the sensitivity is dependent on exposure (see text). 

In Figure |4] we show the resulting SNIa contribution to the 
CGB, for the deflagration model W7 (solid) and the delayed- 
detonation model 5p0z22.23 (dot-dashed). For W7, the shad- 
ing show the uncertainty due to SNIa rate. We see that 
the SNIa contributions are at least a factor ^ 5 smaller 
than current CGB measurements, and up to ^20 depend- 
ing on the SNIa rate and SNIa model. Although the spec- 
trum per SNIa has line features, these are washed out due 
to redshift. Our results are comparable to those of previ- 
ous studies. Altho ugh early studies showed large contri- 
butions from SNIa (Clavton & Silk' '1969; 'Clavton & Ward' 
1975; The et al. 1993; Zdz iarski 1996; Watanabe et al. 1999; 
Ruiz-Lapuente et al.l 200 ll). la ter studies report contribu tions 
to be ^ 10 ("Stri gari et al.l2005l) and >10 (lAhn et al.l2005b less 
than the measured intensity. 

If the dominant fraction of the CGB could be subtracted 
by future detectors, the SNIa contribution could be de- 
tected. The feasibility depends on the true nature of the 
currently measured CGB. Pr oposed sources inc lude various 
populations of AGNs (e.g. , Aiello et al] l2009l) . hot coro- 



nae of AG Ns (jlnoue et aL 2008h, and exotic dark mat- 



ter models (Ahn & Komatsu '2005| ICembranos et all l2007l 
iLawson & Zhitnitskv 2008). The measured CGB may also 
be dominated by detector backgrounds. The true nature is un- 
known, and remains an important quest for future experiments 
to elucidate. Angular-correlatio n techniques may help d iffer- 
entiate the various possibilities (IZhang & Beacomll2004l) . 

3.3. Local SNIa rate measurements 

Knowing the local SNIa rate is a critical prerequisite for 
assessing the prospects for detecting gamma rays from indi- 
vidual SNIa. Many stu dies have discussed th e SNIa rate, most 
no tably in flie 1980s bviOehrels et al.l (ll987h and in the 1990s 
bv lTimmes & Wooslevl(ll997b . which provided much needed 
guidance on gamma-ray detection prospects. Now we have 




847 keV flux [cm"^ s"'^ 

Fig. 6. — The total (i.e. cumulative) number of gamma-ray detectable SNIa 
rate per 10 years, as a function of 847 keV line flux (bottom axis), or distance 
(top axis; assuming 847 keV peak emission of 4.3 X lO*'^ s~^ per SNIa). 
The extrapolation of the cosmic SNIa rate, as well as the rates derived from 
supernova catalogs, are shown and labeled accordingly. The dashed lines are 
3(T detector line sensitivities (10^ s) as labeled. The inset shows part of the 
figure magnified for clarity. 

the advantage of systematic supernova surveys and greatly 
improved SNIa statistics. In addition to supernova catalogs, 
we use our cosmic SNIa rate to arrive at a consistent picture 
for the local SNIa rate. 

More than 5000 supernovae discovered up to the end of 
2009 are listed in the Sternberg As tronomical Institute Su- 
pernova Catalog (SAI; iBartunovet al. 2007), the Asiago Su- 
pernova Catalog, and the catalog maintained by the Central 
Bureau for Astronomical Telegrams. In some cases, the cat- 
alogs disagree on details, but these discrepancies are increas- 
ingly rare for more recent supernovae. In the presence of a 
disagreement, we chose the classification in SAI; only small 
quantitative, and no qualitative, differences appear if prefer- 
ence is given to the other catalogs. We have used the reported 
recession velocities in the SAI catalog for their distances, with 
cross checks for the near est SNIa with catalogs of galaxies 
('Karachent sev et al.ll2004l) . 

SNIa discovered over the most recent 40 years (1970-2009) 
are shown in Figure |5] as a function of distance. We see that 
while rare, there have nonetheless been a number of SNIa 
within 10 Mpc, at approximately 1 per decade (unfortunately, 
none in the recent decade). It is clear that one needs to go only 
a small factor in distance to observe > 10 SNIa per decade. At 
larger distance, we see that the commencement of dedicated 
SNIa searches in the 1990s dramatically increased the number 
of SNIa discovered. Yet these are still underestimates, due to 
missing coverage at large distances. Furthermore the An sky is 
not evenly sampled; the Northern hemisphere is more closely 
observed than the South, resulting in a SNIa discovery ratio 
of approximately 1.4 : 1 (2000-2009, within 100 Mpc). 

The comparison to the cosmic SNIa rate is also revealing. 
The red lines show the extrapolated cosmic SNIa rate, and 
the incompleteness of the catalog is apparent by ^ 30 Mpc. 
With ne xt-genera tion surveys such as the Palomar Transient 
Factory dLaw et al.,.2009.) . SNIa measurements are becoming 
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more complete. Many more SNIa are also being discovered 
pre-maximum, offering more targets for future gamma-ray 
detectors. 

Around 20 Mpc the catalog shows more SNIa than the cos- 
mic extrapolation. The excess is as high as a factor ^ 3 for 
the 20-22 Mpc bin, although this is not statistically strong. 
Even considering the uncertainty in the cosmic SNIa rate this 
excess persists. There is also some excess from the Virgo 
galaxy cluster, located at ~ 18 Mpc in the Northern hemi- 
sphere an d containing in exce ss of 1000 galaxies. The SNIa 
sample of ISmartt et al.l (l2009h . which contains SNIa within 
28 Mpc, similarly shows a high SNIa rate. 

3.4. Review of SNIa gamma-my observations 

There are strong upper limits on the gamma-ray emission 
from individual SNIa. Among the earliest was SN 1986G, 
observed by the Solar Maximum Mission iSMM) instrument 

f[atz& Share* ^ 1990). followed by SN 199 IT (iLichti et alj 
94; Leisine e t"al][l995l:lMorris et al.|[T997l) and SN 1998bu 
sising et al.lll999l) both observed by the Compton Gamma 
Ray Observatory (CGRO). The derived limits depend on the 
fact that SN I986G was subluminous, SN 199 IT was super- 
luminous, and SN 1998bu was normal. Furthermore, dis- 
tance uncertainties weaken the limits. In all cases, limits are 
compatible with t heoretical prediction s within uncertainties, 
though very close dHoflichet al.lll994 . 

On the top x-axis of Figure |5] the distance scale is con- 
verted to a 847 keV number flux, assuming that all SNIa yield 
= 4.3 X 10"*^ s^^ during the peak 10^ s. This is a con- 
servative estimate for a normal SNIa; subluminous and super- 
luminous SNIa would have different axis scalings (Table [T]i. 
The in-flight 3(t narrow-line sensitivity of SPI onboard the 
INTEGRAL satellit e is 3 x 10"^ cm"^ s^^ at - 1 MeV for 
a 10^ s exposure dRoques et alj l2003l) . The equivalent for 
COMPTEL onboard the C GRO satellite is 6 x IQ-^ cmr"^ s^^ 
jSchoenfelder et ailll993h . The Nuclear Compton Telescope 
(NCT) balloon experime nt has a nar row-line sensitivity com- 
parable to INTEGRAL ( Bellm et al.l2009.) the DUAL gamma - 
ray mission is 30 times more sensitive jBoggs et alj|2010h . 
and the proposed Advanced Compton Telescope (ACT) satel- 
lite is 60 times better (Boggs 2006). 

Nearby SNIa are labeled in Figure |5] from which we see 
that those observed by CGRO were just out of range. Note that 
SN 199 IT was a superluminous SNIa, for which the CGRO 
horizon is larger than shown in Figure |5] but still insufficient. 
SN 2003gs was briefl y observed by the INTE GRAL sateflite, 
with no line detection dLeising & Diehlll2009l) . 

The INTEGRAL sensitivity horizon to 847 keV gamma rays 
is '-^ 1 1 Mpc, corresponding to 0. 16 SNIa per year (Figure|6]l. 
The 812 keV flux is smaller, and is visible out to a maximum 
of 4-6 Mpc, corresponding to 0.006-0.03 SNIa per year. The 
847 keV horizon for superluminous and subluminous SNIa 
are ^13 Mpc and ^1 Mpc, respectively. In practice, the en- 
ergy resolution of SPI is better than the Doppler broadened 
width of the lines, so that the signal is spread over multi- 
ple energy bands. This reduces the effectiveness of detec- 
tor background rejection with SPI, and the sensitivity horizon 
decreases by a factor of at most ~ 2 (.Gomez-Gomar et alj 



Il998h . This issue similarly applies to next-generation de- 
tectors which we discuss in the next section. The rates are 
summarized in Table |2] neglecting the reduced detector back- 
ground rejection due to line width. These are also labeled in 
Figure |6] where the cumulative number of catalog SNIa are 
shown as functions of the 847 keV flux (bottom x-axis) and 
distance (top x-axis) 

3.5. Prospects for future SNIa gamma-ray detection 

The initial goal of SNIa gamma-ray studies would be a sin- 
gle long-exposure detection of the 847 keV line from a nearby 
SNIa, which, together with the optical data, would provide a 
robust handle on the amount of ^''Ni synthesized in the explo- 
sion. The prospects for detecting the 847 keV emission from 
^ 1 SNIa per year is not too far from current INTEGRAL sen- 
sitivity (Table 121). 

However, more important in the future is improving beyond 
a single 847 keV detection and measuring the gamma-ray 
light curve, as this ultimately holds the power to distinguish 
between SNIa models. In particular, the period during which 
the system transitions from an optically thick to a thin one 
gives the pcsc (t) in Eq. |5] and provides the best opportunity 
for model discrimination. For example. He-detonation pro- 
duces '^^Ni nearer to the WD surface compared to other mod- 
els, resulting in an earlier transition and hence an earlier rise 
of gamma-ray lines. Deflagration is the other extreme, with 
ignition occurring in the central regions of the WD, resulting 
in slower rise of the gamma-ray emission. These points are 
clear in Figure|2l where the He-detonation models HECD and 
HED8 are shown together with delayed-detonation (W7DT 
and DD202C) and deflagration (W7) models. In addition, 
detection of the 812 keV line from '^^Ni is important, given 
its strong emission during the transition period. In Figure |2l 
both super-luminous SNIa (solid) and normal SNIa (dashed) 
are shown together, but optical observations would distinguish 
between these two classes so we focus on model testing within 
each class. 

First, we show on the right y-axis of Figure [T] the signifi- 
cance that can be achieved by ACT observing a SNIa at 20 
Mpc in 1-day (^ 10^ s) viewing periods. Since MeV gamma- 
ray satellites are dominated by detector backgrounds, the sig- 



TABLE 2 

Per- YEAR SNIa gamma-ray detection rates 



Search Mode 


Sensitivity [cm s 






3 X 10-5 


1 X 10-5 


5 X lO-'' 


847 keV detection [yr"^] 


0.1-0.2 


1.0-1.4 


60-100 


847 keV light curve [yr^^] 


0.1-0.2 


1.0-1.3 


60-80 


812 keV detection [yr^^] 


0.006-0.03 


0.03-0.2 


2-20 


812 keV light curve [yr^^] 


0.001-0.006 


0.005-0.03 


0.5^ 



Note. — Rates are deiived using the cosmic SNIa rate and rates from the 
SAI catalog, whichever is larger The sensitivities 3 X IQ- ^ cm-^ s— ^ and 
5 X IQ-^ <^™~^ correspond to the Scr sensitivity of I NTEGRAL-SP I 
jRoque s et alll2003i) and the baseline Scr sensitivity of ACT )Boggsll2003) . 
respectively. The intermediate sensitivity of 1 X IQ- ^ cm— ^ s— ^ is a hy- 
pothetical satellite for illustrative purposes. The terms "detection" and "light 
curve" are defined as a single lO" s exposure detection and three independent 
lO*" s exposure detections, respectively. The difference between detection 
and light curve prospects reflects the fact that the 847 keV is flatter and 812 
keV is more peaked at their respective peak fluxes. 
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^INTEGRAL SNIa at 20 Mpc - 




Time [days] 

Fig. 7. — Top: 812 keV and 847 keV light curves, for a selection of SNIa 
simulations (see Table[T) at a distance of 20 Mpc. The 3(t detector sensitivi- 
ties (10^ s) are labeled by arrows. On the right axis the statistical significance 
for ACr, over a 1 day viewing period, is shown (see text for how significance 
accumulates with the square-root of the viewing period). ACT will produce 
detailed 812 keV and 847 keV fight curves of all SNIa at 20 Mpc (sublu- 
minous SNIa have not been plotted for clarity, but remain above the ACT 
sensitivity at 20 Mpc). Bottom: the difference between models, taking W7 as 
the reference, shown as the difference in signal count rate per day by ACT. 
The horizontal shading indicates the squai'e root of the background rate. 

nificance of the SNIa signal scales as s = A^sig/A'^bkg^^^ oc 
r^/^, where iVgig is the signal counts, A^bkg is the background 
counts, and T is the viewing period. The targeted Scr sensitiv- 
ity of ACT is 5 X 10~''cm~^ (for 10^ s), and we scale it 
to a viewing period of 1 day. We see that ACT will allow de- 
tailed reconstruction of the 812 keV and 847 keV light curves 
with very large significance at this distance. 

Second, provided that the SNIa is close enough for suffi- 
cient photon statistics, studying the shape of the light curve 
does not depend on knowing the exact distance to the SNIa. 
The shape is distinctly different between models, as we illus- 
trate more clearly in the bottom panel of Figure |2l where we 
show the difference counts at ACT, for a viewing period of 
1 day, all relative to the W7 model. An uncertain distance 
would scale the fluxes accordingly, but retain the shape differ- 
ences between models. For illustration, the square root of the 
background counts is shown as the horizontal shading. ACT 
would easily distinguish between normal SNIa models with 
high significance at this distance. The difference between su- 
perluminous SNIa is smaller than the square root of the back- 
ground counts; better testing would be possible with a coarser 
time binning. 

The physics potential can be scaled to a detector with 
generic properties. The significance shown in Figure [T] (top 
right y-axis) scales as 



^ / 5 X 10-7cm-^s-i \ /lOSy/' 

V (?!>scns J \ Tq J 

where (j)sig is the signal flux, T is the viewing period, and 
0scns is the 3(7 line sensitivity in a To = 10^ s viewing pe- 
riod. Note significance accumulates with the square-root of 
the time period, i.e., 2 consecutive lOcr detection implies an 
overall 10\/2a detection. The square root of the background 
counts (bottom of Figure|7]l scales as 

f 0scns \ / Tq 

where A^s is the effective area of the detector. 

The rate of SNIa within 20 Mpc is ~ 1 per year, and there- 
fore, ACT will strongly constrain SNIa models at a rate of at 
least 1 per year. SNIa further away also yield information, 
with statistical significance gained by coarser time binning. 
At 2-20 times per year, ACT would detect the 812 keV emis- 
sion and build a detailed 847 keV light curve. At 60-80 times 
per year, ACT would detect the 847 keV light curve near peak, 
and at almost 100 times per year, it would detect the 847 keV 
emission and measure the variation in the ^^Ni yield of SNIa. 
We summarize these gamma-ray detection prospects in Ta- 
ble|2l The wide-field all-sky nature of ACr could in principle 
detect SNIa independently of optical discoveries, although in 
practice next-generation optical surveys would bring out the 
full potential of ACT by providing more targets for joint anal- 
ysis. DUAUs focusing capabilities would rely on knowledge 
of the SNIa location from an optical trigger. 

For all the above estimates we have assumed the narrow- 
line sensitivities of detectors as documented. Since SNIa 
gamma-ray lines are expected to be Doppler-broadened by up 
to 3%, this treatment is optimistic, although it does allow us to 
compare the maximum performance of detectors. In practice, 
the broad-line sensitivities are only marginally worse than the 
narrow-line sensitivities. For example, the 3% broaden ed line 
sensiti vity of ACT is 1.2 x 10"^ cm^^ g-i jjj iqG g. Boggsl 
(l2006h ). a factor 2 different from the narrow-hne sensitivity. 
The gamma-ray horizon therefore decrease by ^ 1/ a/2, and 
the numbers in the final column of Table|2]decrease to 20-30, 
16-21, 1-5, and 0.1-1, respectively. Although smaller, these 
rates would still allow rapid progress in the understanding of 
SNIa. Similarly, the significance and y^N^kg for the broad- 
line case can be calculated to be « 7.3 and « 120 in a 1-day 
bin, respectively. The significance remains high, and ^jNhkg 
remains smaller than the model differences, demonstrating the 
physics potential. 

3.6. Discussions 

We focused on the 812 keV and 847 keV lines because they 
are expected to have the highest fluxes. There are other lines 
for which there are interesting prospects, for example the 158 
keV line from ^^Ni and the 1238 keV line from ^^Co. The 
detectability of these lines is only marginally less than that 
of the 812 keV and 847 keV Unes (Eqs.[3}S- Additionally, 
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the 511 ke V line from annihilation of positrons produced in 
^^Co decay (Eq. |4|i is particularly important for understand- 
ing the positron escape fraction, thought to be ~ 1%, but 
quite uncertam jMilne et alJboOll; iLair et alj|2006l) . It has 
important implicati ons for the origin of the G alactic 511 ke V 
emission (see, e.g.. i Beacom & Yukselll2066l and references 
therein). It should also be mentioned that other radioactive 
nuclei such as '^"'Ti and ^^Ni are expected in SNla. The half 
life of ''^Ti ''^Sc is 68 years, implying historical SNla 
as prime targets rather than the SNla discussed in this paper. 
The ^''Ni ^"^Co ^""Fe decay chain, with half lives of 52 
hours and 391 days, would be a more suitable target. How- 
ever, as the production of ^^Ni is not as abundant as ^^Ni in 
most cases, it would be a useful line only for exceptionally 
close SNla. 

Other model-distinguishing features that we have not men- 
tioned include line shift and line width evolutions. During the 
first few weeks, the high optical depth means that much of the 
gamma-ray emission originates from the approaching ejecta, 
with photons escaping essentially radially, so that the lines are 
blue-shifted. Together with the line width, these are measures 
of the expansion velocity, and their time evolution tests SNla 
models. The caveat with the former is that shifts occur mostly 
during the optically thick regime, so model testing is diffi- 
cult except for exceptional cases (e.g., the DD scenario we 
mention next). The line width offers a more promising test, 
in particular with detectors with spectral resolution of order 
\/SX - 200-300 (see, e.g.. iHoflich etaDll998l: iMilne et alJ 
|2004 for detailed predictions). 

The gamma-ray emission for the DD scenario experiences 
significant suppression due to the envelope resulting from the 
merging WDs. In the det2env4 model of Hoflich et al. 
the enveloping matter means that the line shifts remain high 
for a long period of time, up to ^ 10 days, whereas SD sce- 
nario models all fall in a few weeks jHoflich et al. | |l998l) . In 
addition, the envelope strongly suppresses the 812 keV from 
^^Ni decay, so that the flux is weaker than even the deflagra- 
tion model W7. The 847 keV emission also peaks at a later 
stage. Recently, several authors have explored the collision 
of two WDs in dense stellar systems as an alt ernate pathway 
to SNla within the fami ly of DD scenarios jRosswog et alj 
l2009HRaskin et al.ll2009l) ." Shock-triggered thermonuclear ex- 
plosions in such events are predicted to yield ^^Ni masses 
that are sufficient to power subluminous SNla. Although de- 
tailed gamma-ray light curves have not yet been published, 
they could be targets for next-generation gamma-ray detec- 
tors if their gamma-ray emissions are comparable to those of 
subluminous SNla. 

3.6.1. SNla Variants 

Growth of the SNla sample has led to variants of SNla be- 
ing discovered. Deep searches have revealed peculiar faint 
explosions: SN 2002bj, observed by the LOSS survey, was a 
faint SNla showing unusually rapid evolution on time scales 
of days (Poznanski et al. 2010). The low luminosity and short 
rise time (< 7 days) translate to O.15-O.25M0 of ^"^Ni. An- 
other supernova, SN 2005E, has an estimated rise time of 7- 
9 days, very small estimated ^^Ni yield, and spectra show- 



ing a bundance of helium burning products dPoznanski et alj 
I2OIOI) . Its location in an isolated galaxy with little star for- 
mation activity suggests an old progenitor, e.g., WD bina- 
ries, instead of the core collapse of a massive star These re- 
semble so-called ".la" explosions, which has been proposed 
to occur in binary WDs undergoing helium mass transfer 
jBildsten et al. 2007). As the binary evolves, the mass re- 
quired for unstable helium burning increases, until a final 
flash that leads to a faint thermonuclear explosion that is 
one-tenth as bright for one-tenth the time of a normal SNl a 
dShen & Bildstenll2009t IKver et alJl2009l: IShen et alJl2O10h . 
The predicted rates match those observed, both being in the 
range of a few percent of SNla. 

Although the ''^Ni yield is lower than in a normal SNla, 
this is compensated by a larger escape probability for gamma 
rays. Indeed, the rapidly falling light curves of SNla variants 
suggest that of order unity of the ^^Co decay gamma rays are 
escaping already by ^ 20 days. This contrasts with normal 
SNla, where up to 90% of the gamma rays at 20 days are 
down-scattered in the ejecta. It turns out the low ^^Ni yield 
and high escape fraction almost compensate each other. The 
time-integrated escape fraction of gamma-ray photons in nor- 
mal SNla is about 50%, so that SNla variants, with half as 
much ^^Ni and up to twice the total escape fraction, provide a 
comparable gamma-ray output to normal SNla. The low rate 
of SNla variants therefore make them subdominant contribu- 
tors to gamma-ray observations. 

The situation is potentially more interesting for positrons. 
Adopting the cano nical escape fraction 1 % for normal SNla 
(lMilneetalJll999h . SNla variants could contribute up to ^ 2 
times more than normal SNla towards the positron budget, 
depending on their rate. This would help narrow the gap be- 
tween the supernovae positron yield and that required from 
observations of the Galactic 511 keV line. Again, line de- 
tection will prove important: late-time detection of the 511 
keV line would provide a direct constraint of the fraction of 
positrons escaping through the ejecta. 

Another SNla variant is the extremely luminous super- 
Chandrasekhar mass SNla. They exhibit higher magnitudes 
and slower ejecta velocities compared to normal SNla, with 
predicted progenitor masses that exceed the Chandrasekhar 
mass. For example, SN 2003fg was a My = —19.9 
SNla with a predicte d Ni yield of ^ 1. 3 Mq and progeni- 
tor mass of ^ 2 Mq jHowell et al.ll2006l) . The recent super- 
Chandrasekhar mass SN 2009dc has a Ni yield of 1.6 Mq 
( Yaman aka et al. 2009), and spectropolarimetry observations 
suggest the explosion was near s pherical, supporting a truly 
super-Chandrasekhar progenitor ("T anaka et al.ll20 1(y). Rapid 
rotation may support such a massive WD. Binary WDs could 
also produce super-Chandrasekhar mass SNla. If the light 
curves of candidate super-Chandrasekhar mass supernovae 
are indeed powered by a larger-than-normal amount of pro- 
duced nickel, this will give a large gamma-ray signal; if they 
are instead powered by circumstellar interactions, this will not 
give such a large gamma-ray signal. 

4. CONCLUSIONS 
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Gamma rays provide unique clues to the currently debated 
progenitor properties and explosion mechanisms of Type la 
supernovae (SNla). They directly probe the power source of 
SNla and provide tomography of the SNla ejecta. As the 
importance of SNla in astrophysics and cosmology contin- 
ues to grow, the detection of gamma rays become increas- 
ingly essential. The detectability relies on there being suf- 
ficient optical SNla discoveries within the "sensitivity hori- 
zon" of gamma-ray telescopes. In this paper we investigate 
the prospects for studying SNla physics using current and fu- 
ture gamma-ray detectors. Below we summarize our results. 

4.1. Results on SNla rates and SNla progenitors 

We first investigate the SNla rate, which is a prerequisite 
for SNla gamma ray detection prospects. It is also interesting 
since the SNla rate with respect to its progenitor formation 
rate depends on what the SNla progenitors are. We jointly 
analyze the cosmic star formation rate, cosmic SNla rate, and 
the SNla rate derived from SNla catalogs. We deduce a delay 
time distribution (DTD) and SNla fraction that fit the data, 
and discuss the local (< 100 Mpc volume) SNla rate. 

• Delay-times: When SNla rate measurements are cat- 
egorized according to sample size and fraction of 
spectroscopically-identified SNla, we find that the more 
reliable measurements collectively show significantly 
slower evolution with redshift than the star formation 
rate. The difference is due to the delay between pro- 
genitor formation and SNla: we find that a DTD of the 
form cx i^" with a = 1.0 ±0.3 provide s a good fit. The 
substantially-prompt bimodal DTD of iMannucci et alJ 
(2006 ) and the narrow Gaussian DTD around 3.4 Gyr 
of lDahlen et all 120081) do not fit the global data as well. 

• SNla efficiency: For our DTD, we find an efficiency of 
making SNla of (5 ± 1) x IO^^'Mq^ Assuming that 
the SNla progenitor mass range is 3-8 M0, this equates 
to a SNla fraction of 2.4 ± 0.5% (for the Salpeter initial 
mass function; the dependence on this choice is weak). 

• Implication for local SNla rate: The local SNla rate is 
much higher than previously thought. Supernova cata- 
log entries between 2000 to 2009 reveal on average 40 
SNla per year within 100 Mpc, a significant increase 



from the ~ 2 in the 1980s ( Gehre 



in the 1990s ( Timmes & Wooslevlll997' ). Even so, dis 



s et a 



] ll987h and 5.5 



coveries are still severely incomplete outside about 30 
Mpc. The expected true rate within 100 Mpc is about 
« 100 SNla per year The SNla rate within 20 Mpc is 
^ 1 per year 

4.2. Results on SNla gamma rays and SNla explosions 

The detection of gamma rays from SNla directly tests the 
^^Ni mass inferred from optical observations, and also pro- 
vides tomography of the SNla interior by measuring the time- 
dependent escape probability. In the previous section we 
quantitatively discussed how SNla discoveries are becoming 
more complete, and highlighted 20 Mpc as the distance for 
annual SNla discovery. Our main results on SNla gamma-ray 



detection prospects and the physics potential of gamma-ray 
detectors are as follows. 

• CGB contribution: The SNla contribution to the cos- 
mic gamma-ray background (CGB) is at most 10-20% 
of the CGB flux published by SMM and COMPTEL, 
confirming previous results. The origin of the ^ MeV 
CGB, and whether the SNla contribution can be identi- 
fied, therefore remain an important task to be clarified 
by future gamma-ray observations. 

• Current local SNla prospects: Local SNla gamma-ray 
detection prospects are better than thought a decade 
ago, principally driven by the vastly increased rate of 
SNla discoveries. Current gamma-ray satellites probe 
SNla in the rai-e regime: INTEGRAL probes SNla 
within ^ 10 Mpc, occurring at a rate of « 0.1 SNla 
per year This is only somewha t larger than previ- 
ous estimates of ~ 03 p er year dCiehrels etal.lll987^ 
Ijimmes & W ooslevlll99i . but we are more confident 
about the normalization. 

• Future local SNla prospects: The distance for annual 
SNla discovery (20 Mpc) is only a factor 2 further 
than current horizons (Figures |5] and |6]l. Future detec- 
tors with a line sensitivity of 1 x 10~^ cm^^ s^^ — only 
a factor 3 better than that of INTEGRAL — will cross 
this threshold. The proposed ACT satellite, with a 60 
times better line sensitivity than that of INTEGRAL, 
will probe SNla out to ^ 90 Mpc, translating to a rate 
of 100 SNla per year Improved supernova surveys are 
discovering more SNla and at earlier times, making this 
possible. 

• Implication for explosions: Nearby SNla will be tar- 
gets for detecting gamma ray lines from both ^^Ni and 
^^Co decays, and for their light curve reconstruction. 
The ACr satellite will give a hugely significant > IOOct 
detection of gamma-ray light curves every year from 
SNla within 20 Mpc. These SNla will allow detailed 
analysis of SNla explosion physics, providing unprece- 
dented understanding of the SNla explosion mecha- 
nism. A more modest detector with a line sensitivity of 
1 X 10^^ cm^^ s^^ would measure the 847 keV light 
curve for tomography once a year, and detect the 812 
keV line ^0.1 per year (Table |2]l. 

Given the importance of SNla from cosmology to nucle- 
osynthesis, the need to understand the mechanisms of SNla 
will only increase with time. Detecting their gamma ray emis- 
sion is the only way to probe their inner physics. Judging from 
the SNla rates, even a modest improvement in gamma-ray line 
sensitivity would reach the distance range with annual detec- 
tions. Proposed next-generation gamma-ray satellites are in 
an excellent position to rapidly offer revolutionary results. 
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Type Ia supernova rate measurements 
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Kuznetsova et al. (2008) 
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Note. — Redshifts or distance noted are representative, not exact. All values have been corrected to h = 0.7 and 
to match our (rj^j, f2A)=(0.3, 0.7) cosmology. We quote rates in units of SNu (SNu = S'A'(100 yr)-^ (10^° L|)-i) 
or yr~^ Mpc"^^ whichever is reported by the authors. The quoted errors are statistical followed by systematic; where 
only one is present, the systematic is the one not available. A^ia is the number of SNIa used in the study, and the 
percentage of A'^i^ that are spectroscopically identified as SNIa are shown. 

" Hardin et al. (2000) uses the EROS SN sample (1997), with Nja = 4. To deiive the volumetric value we use 
their SNu value with the 2dF luminosity density, with the error on the luminosity density added in quadrature to the 
systematics. 

" Dilday et al. (2010) uses SDSS-II supernova survey data with a total A^/^, = 516. 



